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This  paper  investigates  the  three-dimensional  turbulent  flow  in  a  chimney  used  for  venting  flue  gas  from 
either  a  boiler  or  stove.  The  thermoelectric  generator  (TEG)  modules  are  embedded  in  the  chimney  walls. 
To  understand  the  power  output  performance  of  the  TEG  module,  three-dimensional  numerical  simu¬ 
lations  combining  convection  and  radiation  effects,  including  the  chimney  tunnel,  TEG  modules,  plate-fin 
heat  sinks  and  cold  plates,  based  on  water  cooling  are  developed  and  solved  simultaneously.  The  effects 
of  operational  parameters  such  as  the  flue  gas  velocity  (Vjn  =  3,  5  and  10  m/s)  and  flue  gas  temperatures 
(Tgas  =  500,  600  and  700  K)  on  the  flow  and  heat  transfer  are  determined.  The  influences  of  the  plate-fin 
height  (Hfin)  and  number  of  fins  (N),  ranging  from  0  mm  <  Hfjn  <100  mm  and  4  <  N  <  8  respectively,  on 
the  power  output  and  pressure  drop  are  also  described  in  detail.  It  is  worthy  of  note  that  the  net  electric 
power  (Pnet)  of  the  TEG  module  was  obtained  using  the  ideal  electric  power  (Pteg)  minus  the  extra 
pumping  power  (Pfan).  The  numerical  results  for  the  power  vs.  current  (P— I)  curve  are  in  good  agreement 
with  the  experimental  data  within  an  error  of  9%. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  effects  include  the  Seebeck  effect,  the  Peltier  ef¬ 
fect  and  the  Thomson  effect,  and  they  are  also  accompanied  by  other 
effects,  such  as  the  Joule  heating  and  Fourier  heat  conduction  laws. 
Seebeck  effects  in  semiconductor  materials  induce  electromotive 
force  due  to  temperature  gradients,  and  Peltier  effects  induce  tem¬ 
perature  gradients  when  a  DC  electrical  current  is  applied.  Addi¬ 
tionally,  thermoelectric  generators  present  potential  applications  in 
the  conversion  of  low  level  thermal  energy  into  electrical  power. 
Especially  in  the  case  of  waste  heat  recovery,  it  is  unnecessary  to 
consider  the  cost  of  thermal  energy  input,  and  there  are  additional 
advantages,  such  as  energy  savings  and  emission  reduction,  so  the 
low  efficiency  problem  is  no  longer  the  most  important  issue  that 
needs  to  be  taken  into  account  [1].  In  general,  a  thermoelectric 
generator  (TEG)  consists  of  a  number  of  semiconductor  pairs  that  are 
connected  electrically  in  a  series  and  thermally  in  parallel,  and  each 
pair  includes  a  p- type  and  an  n-type  element.  Although  in  theory, 
a  single  piece  of  semiconductor  material  could  work,  a  series 
connection  is  used  to  meet  the  high  voltage  potential  requirements. 
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P-type  and  n-type  elements  are  alternated  to  assure  that  the  carriers 
transport  in  the  same  direction. 

The  reported  performance  data  is  usually  analyzed  using  con¬ 
ventional  non-equilibrium  thermodynamics  [2,3].  With  reference 
to  the  structure  of  a  TEG,  it  has  been  demonstrated  that  a  significant 
increase  in  the  power  output  from  a  module  can  be  achieved  by 
modifying  the  geometry  of  the  thermoelectric  elements  [4,5].  In 
subsequent  research,  Rowe  [6,7]  provided  efficiency  in  a  couple  of 
solar-powered  TEGs  and  reviewed  applications  of  nuclear-powered 
TEGs  in  space.  Chen  et  al.  [8]  used  an  irreversible  model  to  study  the 
performance  of  a  TEG  with  external  and  internal  irreversibility.  The 
optimal  range  of  the  parameters  for  the  device-design  was  deter¬ 
mined,  and  the  problems  relative  to  the  maximum  power  output 
and  maximum  efficiency  were  discussed.  Xuan  et  al.  [9]  employed  a 
phenomenological  model  to  study  the  effects  of  internal  and 
external  interface  layers  on  thermoelectric  module  performance. 
Sets  of  general  performance  formulas  have  also  been  derived. 
Pramanick  and  Das  [10]  performed  a  study  on  the  structural  design 
of  a  thermoelectric  module.  A  model  for  a  cascaded  TEG  that  is 
based  on  finite  time  thermodynamics  has  also  been  developed. 
Astrain  [11]  used  the  numerical  method  of  finite  differences  to 
evaluate  the  performance  of  the  TEG  system,  including  the  heat 
exchangers,  the  heat  source  and  the  heat  sink. 
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Nomenclature 

Greek  symbols 

a 

Seebeck  coefficient  (V  K_1) 

Ab 

heat  sink  base  area  (m2)  A  =  Lfin  x  Lfin 

p 

Thomson  coefficient  (V  K"1) 

E 

electric  field  (V  m-1) 

e 

dissipation  rate 

h 

effective  heat  transfer  coefficient  (W  m-2  K-1) 

ee 

extinction  coefficient  (m-1) 

Hfin 

fin  height  (mm) 

V 

TEG  conversion  efficiency  (%)  p  =  Pteg/Q/i 

I 

electric  current  (A) 

K 

turbulent  kinetic  energy 

J 

electric  current  density  (A  m-2) 

P 

electrical  resistivity  (Q  m) 

k 

thermal  conductivity  (W  m-1  K_1) 

O 

Stefan  Boltzmann  constant  (W  m-2  K  4) 

fin  length  (mm) 

(7a 

absorption  coefficient  (m-1) 

n 

pairs  of  p- type  and  n-type  semiconductor 

(7s 

scattering  coefficient  (m_1) 

p 

power  output  (W) 

0J 

uncertainty 

P/Ab 

power  density  (W  m-2) 

K 

general  dependent 

Pr 

Prandtl  number 

Q 

heat  transfer  rate  (W) 

Subscripts 

<7 

heat  generation  (W  m-2) 

c 

cold  side 

R 

residual  sum 

cer 

ceramic  substrate 

Rl 

external  load  resistance  (Q) 

cond 

conductive  copper 

ReL 

Reynolds  number  based  on  fin  length 

/ 

fluid 

5fin 

fin  spacing  (mm) 

fan 

fan 

Steg 

TEG  module  spacing  (mm) 

gas 

flue  gas 

T 

temperature  (I<) 

h 

hot  side 

tb 

base  thickness  (mm) 

net 

net 

ffin 

fin  thickness  (mm) 

P-N 

p-  and  n-type  thermoelectric  semiconductor 

V 

electric  potential  (V) 

s 

solid 

Vin 

gas  velocity  (m  s-1) 

TEG 

thermoelectric  generator 

vmp 

voltage  at  maximum  power  (V) 

water 

cooling  water 

Voc 

open  circuit  voltage  (V) 

x,  y,  z 

Cartesian  coordinates 

In  regard  to  the  application  of  thermoelectric  technology  in 
industry,  there  has  been  a  growing  interest  in  TEGs  using  various 
heat  sources  in  power  plants  for  combustion  of  waste,  geothermal 
energy,  and  other  industrial  heat-generating  processes  [12,13]. 
Crane  and  Jackson  [14]  studied  numerical  heat  exchanger  models 
integrated  with  models  for  Bi2Te3  thermoelectric  modules  which 
were  validated  against  experimental  data  from  previous  cross-flow 
heat  exchanger  studies  as  well  as  from  experiments  using  ther¬ 
moelectric  modules  between  counter-flow  hot  water  and  cooling 
air  flow  channels.  Dai  et  al.  [15]  studied  the  possibility  of  a  TEG 
system  based  on  liquid  metal  which  served  to  harvest  and  transport 
waste  heat.  The  experimental  results  for  the  TEG  system  were 
discussed,  and  a  calculated  efficiency  of  2%  in  the  entire  TEG  system 
was  obtained.  Furthermore,  in  the  application  of  the  waste  heat 
recovery  power  generation,  there  have  been  many  conceptual  de¬ 
signs  for  a  power  conversion  system.  Thacher  et  al.  [16]  investi¬ 
gated  the  feasibility  of  waste  heat  recovery  from  exhaust  in  a  light 
truck  by  connecting  a  series  of  16  TEG  modules,  which  showed 
good  performance  at  high  speeds.  Hsiao  et  al.  [17]  constructed  a 
mathematical  model  to  predict  the  performance  of  a  TEG  module 
attached  to  a  waste  heat  recovery  system.  It  was  observed  that  the 
TEG  module  presented  better  performance  on  an  exhaust  pipe  as 
opposed  to  a  radiator.  Niu  et  al.  [18]  constructed  an  experimental 
TEG  unit  that  was  used  to  examine  the  influences  of  the  main 
operating  conditions,  hot  and  cold  fluid  inlet  temperatures,  flow 
rates  and  load  resistance  on  both  power  output  and  conversion 
efficiency.  Champier  et  al.  [19]  studied  the  influence  of  the  pressure 
on  the  contact  resistance  between  the  TEG  and  stoves,  and  the 
experimental  measurements  are  compared  with  the  results  of  a 
theoretical  analysis  used  one  dimension  thermoelectric  model. 

Because  the  power  output  of  a  TEG  module  without  any  heat 
collector  is  lower,  to  obtain  high  quality  power  output,  it  is  feasible 
to  insert  a  plate-fin  heat  sink  into  a  chimney  tunnel,  which  can 


increase  the  convective  heat  transfer  area  and  reduce  thermal 
resistance  from  the  flue  gas  to  the  TEG  module.  However,  it  pro¬ 
duces  a  disadvantage  in  that  the  pressure  drop  of  the  system  is 
relatively  increased.  Jonsson  and  Moshfegh  [20]  studied  the  ther¬ 
mal  and  hydraulic  performance  of  42  different  heat  sinks,  including 
plate-fin,  strip-fin,  and  pin-fin  heat  sinks,  by  varying  the  fin  height, 
fin-to-fin  distance,  and  the  tip  and  lateral  clearance.  In  addition,  an 
empirical  bypass  correlation  of  pressure  drop  was  proposed  in  their 
work.  Lee  [21]  developed  an  analytical  simulation  model  for  pre¬ 
dicting  the  thermal  and  pressure  characteristics  of  fin  heat  sinks  in 
a  partially  confined  configuration.  Barrett  and  Obinelo  [22]  studied 
straight-channel  longitudinal  fin  heat  sinks  under  forced  air  cool¬ 
ing.  In  this  report,  the  influence  of  flow  bypass  phenomena  on  the 
thermal  resistance  and  pressure  drop  across  the  heat  sinks  were 
also  considered  and  discussed. 

According  to  the  literature  survey  on  this  topic  [23],  it  is 
recognized  that  most  of  the  previous  theoretical  models  have  been 
limited  to  one  dimensional  problems.  However,  simplified  models 
in  which  the  p—n  element  pair  is  simply  treated  as  a  single  bulk 
material  are  not  sufficient  to  accurately  evaluate  the  detailed 
behavior  of  most  practical  fluid  systems  with  multi-dimensional 
construction,  irregular  geometry,  dynamic  variations  in  flow  field, 
and  complicated  heat  transfer,  nor  can  they  precisely  indicate 
temperature  distributions  and  nonuniform  heat  flux  to  the  TE 
model  as  the  boundary  condition.  Only  a  few  studies  have  been 
conducted  in  which  a  numerical  analysis  of  a  three-dimensional 
TEG  model  has  been  conducted.  Sandoz-Rosado  and  Stevens  [24] 
presented  a  finite  element  model  to  address  three-dimensional 
transport  in  a  TEG  model.  Simultaneously,  a  new  method  of  pre¬ 
dictive  error  mapping  for  one-dimensional  analytical  models  was 
developed  and  demonstrated.  Chen  et  al.  [25]  proposed  a  three- 
dimensional  TEG  model  which  was  implemented  by  computa¬ 
tional  fluid  dynamics  commerce  software  (FLUENT)  and  its 
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Fig.  1.  Schematic  of  the  thermoelectric  module  with  a  plate-fin  heat  sink  for  waste 
heat  recovery. 


user-defined  functions.  Cheng  et  al.  [26]  provided  a  three- 
dimensional  numerical  model  to  predict  the  transient  thermal 
behavior  of  TECs.  The  results  showed  that  both  the  temperatures  of 
the  hot  and  the  cold  ends  increase  with  the  cooling  load  and  that 
the  value  of  the  COP  (coefficient  of  performance)  linearly  increases 
with  the  cooling  load.  Jang  et  al.  [27]  used  a  numerical  optimization 
technique  in  a  geometrical  optimization  for  a  TEG  module  spreader 
to  obtain  efficient  power  output  balancing  with  reasonable  mate¬ 
rial  costs.  The  results  showed  the  optimum  spreader  length 
and  thickness  with  the  fixed  module  spacing  for  Pareto-optimal 


configurations.  Jang  and  Tsai  [28]  studied  the  optimization  of  TEG 
module  spacing  and  its  spreader  thickness  using  a  simplified 
conjugate-gradient  method.  In  this  report,  the  effects  of  tempera¬ 
ture  difference  and  waste  gas  heat  transfer  coefficients  were  dis¬ 
cussed,  and  the  numerical  data  were  in  good  agreement  with  the 
experimental  data. 

Based  on  the  above  literature  review,  up  to  now,  there  have  been 
no  papers  concerning  the  three-dimensional  turbulent  flow  and 
heat  transfer  in  a  waste  heat  recovery  system  with  a  TEG  module 
coupled  with  a  plate-fin  heat  sink,  as  shown  in  Fig.  1.  This  fact  has 
motivated  the  present  study.  The  purpose  of  the  present  study  is  to 
build  up  a  complete  computational  domain,  including  the  chimney 
tunnel  flow,  the  TEG  module,  and  the  plate-fin  heat  sink,  and  to 
directly  solve  the  three-dimensional  turbulent  flow  for  an  integral 
waste  heat  recovery  system.  Because  of  the  high  temperature  flue 
gas  flowing  into  the  chimney  tunnel,  the  radiation  effect  is 
considered.  Numerical  simulations  focus  on  the  effects  of  the 
geometric  parameters  of  the  plate-fin  (fin  height,  Hfin  and  number 
of  fins,  N)  and  the  operational  parameters  (frontal  velocity,  Vin  and 
waste  gas  temperatures,  Tgas)  on  the  TEG  module  power  output  and 
on  the  chimney  tunnel-induced  pressure  drop.  Furthermore,  an 
experiment  is  executed  to  verify  the  accuracy  of  the  numerical 
results  for  the  TEG  power  output. 

2.  Mathematical  analysis 

A  schematic  diagram  of  thermoelectric  modules  with  a  plate-fin 
heat  sink  is  shown  in  Fig.  2.  A  thermoelectric  module  consists  of  a 


Fig.  2.  The  computational  domain  of  the  thermoelectric  module  (a)  dimensions  of  the  thermoelectric  module  (b)  variables  of  the  plate-fin  heat  sinks. 
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Table  1 

Geometrical  parameters  and  thermal— electrical  properties  of  the  TEG  module  [29]. 


Geometrical  and 
thermal-electrical  properties 

Symbol 

Value 

1.  Plate-fin  heat  sink 

Thermal  conductivity 

ksp 

155  (aluminum) 

(W  nr1  K1) 

Fin  thickness  (mm) 

ffin 

1.5 

Fin  length  (mm) 

f-fin 

30 

Fin  base  thickness  (mm) 

tb 

5 

2.  Ceramic  substrate 

Thickness  (mm) 

Hcer 

0.75 

Thermal  conductance 

kCer 

49.2 

(W  m1  K-1) 

3.  Conductive  copper 

Thickness  (mm) 

ffcond 

0.12 

Thermal  conductance 

kcond 

403 

(W  ITT1  K-1) 

Seebeck  coefficient  (V  IC1) 

^cond 

1.4  x  10"5 

Resistivity  ( Q  m) 

Pcond 

1.44  x  10"8 

4.  P— N  thermoelectric 
elements  (for  one  pair) 

Thickness  (mm) 

Hp-n 

1.68 

Thermal  conductance 

kp-N 

1.6 

(W  m1  K-1) 

Seebeck  coefficient  (V  IC1) 

ap-N 

±2.0  x  10"4 

Electrical  resistivity  (Q  m) 

Pp-n 

9.0  x  10“6 

plate-fin  heat  sink,  a  TEG  and  a  cold  plate  based  on  water  cooling. 
The  dimensions  of  the  assembled  square  TEG  unit  are 
30  mm  x  30  mm  x  3.42  mm,  and  other  detailed  geometrical  pa¬ 
rameters  and  thermal-electrical  properties  of  the  TEG  module  are 
summarized  in  Table  1  [29].  The  major  components  of  flue  gas  from 
a  steel  plant  are  N2  (78%),  C02  (13%),  H20  (6%),  02  (2.6%),  CO  (0.2%) 
and  S02  (0.2%).  The  thermophysical  properties  of  the  flue  gas  can  be 
calculated  as  in  Irvine  and  Liley  [30],  as  shown  in  Table  2.  A  typical 
thermoelectric  generator  consists  of  many  positive-type  (p- type) 
and  negative-type  (n-type)  semiconductor  legs  that  are  connected 
by  copper  conducting  strips  in  serials.  However,  a  basic  thermo¬ 
electric  element  is  composed  of  both  a  p- type  and  an  n-type 
semiconductor  leg.  Since  we  use  the  thermoelectric  models  on  the 
basis  of  an  equivalent  current  method,  which  means  there  is  an 
equal  electric  current  and  an  equal  cross  area  of  thermoelectric 
elements,  only  one  pair  is  considered  here.  In  the  present  study,  a 
TEG  module  with  31  pairs  of  p- type  and  n-type  semiconductors  is 
simplified  to  one  pair.  Then,  the  equivalent  thermal-electrical 
properties  including  thermal  conductivity  /<p_n,  electrical  resistivity 
Pp-n  and  Seebeck  coefficient  ap-N  are  listed  in  Table  3. 


Table  2 

The  thermophysical  properties  of  the  flue  gas  [30]. 


Thermophysical  property 

Tgas  =  500  K 

Tgas  =  700  K 

Thermal  conductivity  (W/m  K) 
Viscosity  (N  s/m2) 

Specific  heat  (J/kg  K) 

Density  (kg/m3) 

3.152  x  10"5 
2.171  x  10"5 
1.0796 

0.9025 

3.804  x  10"5 
2.553  x  10"5 
1.1039 

0.7261 

Table  3 

The  equivalent  thermal— electrical  properties  when  the  TEG  module  with  n  pairs  of 
p- type  and  n-type  semiconductors  is  simplified  to  one  pair. 

Thermo-electrical  property 

p- type  and  n-type 
(for  one  pair) 

When  TEG  with  n 
pairs  is  simplified 
to  1  pair 

Thermal  conductivity  (W/m  K) 
Electrical  resistivity  (Q  m) 
Seebeck  coefficient  (V/I<) 

kp-N 

Pp-n 

ap-N 

kP_N (unchanged) 

n2pp~N 
n  ap-N 

2  A.  Governing  equations  of  thermoelectricity 

In  a  thermoelectric  analysis,  the  three-dimensional  governing 
equations  of  heat  and  electric  current  in  TEG  at  steady  state  can  be 
cast  into  the  following  general  form: 


V-<7  =  <7  (1) 

V-T  =  0  (2) 

In  the  foregoing  equation,  the  heat  flux  vector  q  in  terms  of  the 
electric  current  density  vector  J  and  Fourier’s  law  fcvT  are 
expressed  as  the  following: 


q  =  aTJ  -  kVT  (3) 

7  =  I(i-«vr),  (4) 

where  E  =  -W  is  the  electric  field  intensity  vector,  and  V  is  the 
electric  scalar  potential.  When  substituting  Eq.  (3)  into  Eq.  (1),  the 
governing  equation  of  heat  flow  becomes 

V  •  (fcVT)  —  V  •  (aT  +  q  =  0.  (5) 

The  heat  generation  termj  shown  in  the  preceding  equation 
includes  the  electric  power  J  ■  E  spent  on  Joule  heating  and  on 
work  against  the  Seebeck  field  ctVT.  In  Eqs.  (4)  and  (5),  /<,  p  and  a  are 
the  thermal  conductivity,  electrical  resistivity  and  Seebeck  coeffi¬ 
cient,  respectively.  Furthermore,  with  the  assumptions  of  constant 
material  properties  k  and  p,  the  temperature  distribution  in  the 
interior  of  TEG  can  be  described  by  a  linear  second-order  differ¬ 
ential  equation  as 

kV2T  -  pj VT  +  pj2  =  0.  (6) 

Physically,  the  terms  appearing  on  the  left-hand-side  of  the 
above  equation  represent  the  heat  conduction,  Thomson  heat  and 
Joule’s  heat,  respectively.  However,  the  Thomson  effect  is  neglected 
because  the  Seebeck  coefficient  is  independent  of  temperature  in 
the  present  study,  and  thus  Eq.  (6)  can  be  expressed  as 
-±2 

kW2T  +  pJ  =0.  (7) 


adiabatic 


Fig.  3.  Boundary  condition  of  a  thermoelectric  generator  module. 
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wall 

(adiabatic) 


TEC  module 


wall 

periodic  (convection) 


wall  (convection) 


for  mass,  momentum  (Reynolds-averaged  Navier-Stokes  equa¬ 
tion),  energy,  turbulent  kinetic  energy,  /c,  and  dissipation  rate,  ey  in 
the  fluid  region  can  be  expressed  in  the  tensor  form  as  follows: 


M  =  0 

dXi 


(8) 


_a 

ax. 


.P^Uj) 


dp_ 

dXi 


dXj 


Meff 


dUi  , 


dXj 


dU{ 

dpq  I  ~~  pU'iUi 


(9) 


Apcp(n,T)  =^+^-^(fcf-PQ,^)-v-?r  do) 


a  a  / jueff  dK  \ 

=-wX-t^i)+p{PT-£) 


(id 


Fig.  4.  Boundary  condition  and  computational  grid  system  of  a  chimney  tunnel  with  Q  _ 
four  TEG  modules.  (Pui£) 


_9_  /^eff  9e\ 

dXi  \  (Te  dxJ 


(12) 


Table  4 

Comparisons  between  our  1  -D  TEG  numerical  solution  and  previous  investigators 
[42]  for  Th  =  500  K  and  Tc  =  300  K. 


Voc  (volt) 

Pmax  (W) 

Present  simulation 

0.023 

0.03174 

Fraisse  et  al.  [42] 

0.0236 

0.0326 

Error  (%) 

2.54 

2.64 

where 


Pr 


k2 

Meff  =  M  +  Mt,  Mt  =  pCii— 


=  0.09,  Ci  =  0.15,  C2  =  1.90,  C  =  0.25,  ak 
=  0.75,  and  ae  =  1.15. 


2.2.  Governing  equations  of  flow 

In  the  present  study,  simulations  are  performed  for  Reynolds 
numbers  Rei  (based  on  the  plate-fin  length  (Lfin)  in  a  flow  direction 
of  30  mm  and  a  frontal  velocity  1-10  m/s)  ranging  from  1000  to 
10,000.  Therefore,  the  fluid  is  considered  incompressible  with 
constant  properties,  and  the  flow  is  assumed  to  be  a  three- 
dimensional,  steady,  turbulent  flow.  The  dimensionless  equations 


Equation  (9)  contains  Reynolds  stresses  that  are  modeled  by 
Chen’s  extended  k— e  turbulence  model  (Chen  [31  ],  Wang  and  Chen 
[32]).  In  Chen’s  model,  both  the  production  time  scale  and  the 
dissipation  time  scale  are  used  in  closing  the  e  equation.  This  extra 
production  time  scale  is  claimed  to  allow  the  energy  transfer 
mechanism  of  turbulence  to  respond  to  the  mean  strain  rate  more 
effectively.  This  results  in  an  extra  constant  in  the  e  equation.  As  to 
the  velocity  distribution  in  the  near- wall  region  (y+  <  11.63),  the 
following  law  of  the  wall  (Liakopoulos  [33],)  is  applied: 


Fig.  5.  Schematic  diagram  of  the  experimental  setup. 
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l(y+2  -7.37y+  +  83.3)°'79_ 

+  5.63tarr1(0.12.y+-0.44l)  -3.81,  (13) 

where, 

y+  =  P}W  and  ut  =  ,/±.  (14) 

P  VP 

Turbulent  kinetic  energy,  /c,  and  dissipation  rate,  e,  in  the  near 
wall  region  can  be  expressed, 

K|y+  =  u?qr1/2  and  £|y+  =  (15) 


2.3.  Governing  equations  of  gas  radiation  modeling 

In  the  energy  equation  Eq.  (10),  in  addition  to  the  conductive 
and  convective  terms,  the  radiative  term  as  the  divergence  of  the 
radiative  heat  flux  is  also  included.  The  radiative  source  term 
V  •  q  r  is  obtained  by  integrating  the  radiation  intensity  fields  over 
all  possible  directions  and  wavelengths: 


4tc 


In  the  above  equation,  oa  is  the  absorption  coefficient  of  the  flue 
gas.  By  a  careful  correlation  of  several  sources  of  experimental  data 
from  Hottel  [34],  the  values  of  the  gas  absorption  coefficient 
(da  =  0.44)  are  estimated.  /(  r ,  s )  is  the  spectral  radiation  intensity 
at  the  situation  r  and  in  the  direction  s  .  To  obtain  this  term,  it  is 
necessary  to  solve  the  radiative  transfer  equation.  This  equation  for 
an  absorbing,  emitting  and  scattering  gray  medium  can  be 
expressed  as  [35]: 

(s-v)/(r-s)  =  -eelp-s)  +  <xa/b(r) 

+  il  (17) 

4tt 

where  as  is  the  scattering  coefficient;  ee  =  da+ds  is  the  extinction 
coefficient,  and  0(  s ,  s  )  is  the  scattering  phase  function  for  the 
radiation  from  incoming  direction  s  confined  within  the  solid 
angle  dQ'  to  scattered  direction  s  confined  within  the  solid  angle 
dQ.  Since  the  scattering  coefficient  of  gas  molecules  is  negligible  in 
comparison  to  the  absorption  coefficient,  we  neglect  the  effect  of 
scattering  in  this  calculation,  that  is,  as  =  0.  In  this  study,  the 
radiative  transfer  equation  is  solved  numerically  by  the  discrete 
ordinates  method  (DOM)  [36,37]. 


Table  5 

Summary  of  estimated  uncertainties. 


Primary  measurements 

Derived  measurements 

Parameter 

Uncertainties 

Parameter 

Uncertainties 

Uncertainties 

Regas  =  1000 

Regas  =  10,000 

Vin 

0.3-1% 

Itegas 

±1.2% 

±0.37% 

JTlwater 

0.5% 

Qgas 

±3.5% 

±1.5% 

I 

0.2% 

hh 

±5.6% 

±2.1% 

V 

0.1% 

P 

±3.2% 

±3.2% 

Tgas 

0.1  °C 

f water 

0.05  °C 

2.4.  Boundary  conditions 

Since  the  governing  equations  are  elliptic  in  regard  to  spatial 
coordinates,  it  is  necessary  to  impose  boundary  conditions  at  all  of 
the  boundaries  in  the  computational  domain.  The  boundary  con¬ 
ditions  and  computational  domains  of  the  TEG  module  and  the 
chimney  tunnel  are  illustrated  in  Figs.  3  and  4,  respectively.  For 
the  TEG  module,  the  constant  heat  transfer  coefficient  hc  is 
specified  at  the  cold  side  of  the  TEG,  and  it  is  estimated  using  the 
Dittus— Boelter  equation  [38].  In  addition,  the  other  remaining 


Fig.  6.  V-l  and  P-I  curves  of  experimental  and  numerical  data  with  Vin  =  5  m/s  and 
Tgas  =  600  K. 
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surfaces  outside  the  chimney  are  set  to  be  adiabatic.  Our  numerical 
analysis  were  carried  for  the  entire  range  of  external  load  resistance 
from  Ri  =  0  (open  circuit)  to  Ri  =  oo  (short  circuit).  For  an  open 
circuit  situation,  the  voltage  at  the  p- type  thermoelectric  element  is 
V  =  0  V,  and  the  current  density  at  the  n-type  thermoelectric 
element  is  J  =  0  A/m2. 

The  upstream  boundary  of  the  computational  domain  is  estab¬ 
lished  at  a  distance  of  100  mm  in  front  of  the  leading  edge  of  the 
first  TEG  module.  Uniform  flow  with  the  velocity  Vin  (1—10  m/s)  and 
temperature  T[n  (500  K,  600  K  and  700  K)  are  assumed.  For  tur¬ 
bulent  flow  calculation,  the  inlet  turbulence  intensity  is  calculated 
from  the  following  equation  and  is  around  3%. 


where  u!  is  the  root-mean-square  of  the  turbulent  velocity  fluctu¬ 
ations,  which  is  obtained  from  the  hot-wire  measurement,  and  uavg 
is  the  mean  velocity. 

At  the  downstream  end  of  the  computational  domain,  located 
seven  times  the  tube  diameter  from  the  last  downstream  row  tube, 
the  streamwise  gradients  (Neumann  boundary  conditions)  for  all 
the  variables  are  set  to  zero.  At  the  solid  interface,  a  no-slip  con¬ 
dition  is  specified.  Due  to  the  symmetry  of  the  flow  domain,  cal¬ 
culations  are  performed  for  a  half  tunnel  width,  and  a  periodic 
match  is  established  on  both  the  top  and  bottom  of  the  domain. 
Additionally,  at  the  solid-fluid  interface  of  heat  sinks,  Ts  =  Tf,  -ks  • 
dTs/dn  =  —kfdTf/dn. 

For  the  radiative  boundary  conditions,  the  walls  are  assumed  to 
emit  and  reflect  diffusely  with  constant  wall  emissivity,  ew  =  0.8 
[39].  In  addition,  the  inlet  and  outlet  sections  are  considered  for 
radiative  transfer  to  be  black  walls  at  the  fluid  temperatures  of  the 
inlet  and  outlet  sections,  respectively. 


2.4.1.  Performance  of  a  TEG 

The  performance  of  a  thermoelectric  generator  can  be  charac¬ 
terized.  The  power  generated  from  a  thermoelectric  generator  can 
be  calculated  as  follows: 

PTEG  =  IV-  (19) 


The  TEG  conversion  efficiency  is  given  by 


where  Qh  is  the  heat  transfer  rate  of  the  hot  side. 


(20) 


3.  Numerical  methods 

Since  there  is  no  analytical  solution  in  this  problem,  the  nu¬ 
merical  simulations  were  carried  out  using  a  commercial  CFD 
software  package  (CFD-RC)  [40].  Fig.  4  illustrates  the  physical 
model  and  computation  grid  system  for  the  TEG  module  with  the 
plate-fin  heat  sink  in  Cartesian  coordinates.  A  grid  system  of  83,722 
grid  points  is  typically  adopted  in  the  computational  domain.  Here, 
the  grid  pattern  is  coarsely  displayed  for  the  convenience  of  the 
reader  to  visualize  the  computational  regions. 

The  finite-volume  method  was  used  to  solve  the  governing 
equations.  It  involves  subdividing  the  region  in  which  the  flow  is  to 
be  solved  into  individual  cells  or  control  volumes  so  that  equations 
can  be  integrated  numerically  on  a  cell-by-cell  basis  to  produce 
discrete  algebraic  (finite- volume)  equations.  Several  algorithms  for 
calculating  the  pressure  field  are  available  in  the  program  (SIMPLE, 
PISO,  SIMPISO).  For  these  steady-state  computations,  the  SIMPLE 
algorithm  [41]  was  found  to  give  the  shortest  computer  running 
times  and  was  used  in  all  of  the  calculations.  Additionally,  the 
relative  error  of  SIMPLE  algorithm  compared  with  the  other 
methods  (PISO,  SIMPISO)  was  within  3%. 

Prior  to  computation,  a  verification  of  the  grid-independence  of 
the  numerical  solution  was  performed  in  order  to  ensure  the  ac¬ 
curacy  and  validity  of  the  numerical  results.  A  grid  system  of  83,722 
grid  points  was  typically  adopted  in  the  computational  domain.  In 
order  to  validate  the  independency  of  the  solution  on  the  grid, 
three  different  grid  systems  were  investigated,  which  included 
about  61,636,  83,722  and  100,954  nodes,  respectively,  for  a  TEG 
module  with  a  plate-fin  heat  sink  in  a  chimney  tunnel.  The  relative 
errors  in  the  local  pressure  (P),  temperature  (T)  and  electric  po¬ 
tential  (V)  between  the  solutions  of  83,722  and  100,954  were  found 
to  be  within  3.0%.  Computations  were  performed  on  an  Intel  Core 


Hfln  =  25  mm  Hfm  =  50  mm  Hfm  =  75  mm 

Temperature 


Electric  potential 
1,41  (V) 


Fig.  7.  The  temperature  and  electric  potential  distributions  of  the  first  thermoelectric  module  with  a  plate-fin  heat  sink  for  N  =  6,  Vin  =  10  m/s  and  Tga s  =  600  K. 
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i7-2600  3.40G  personal  computer,  and  typical  CPU  times  were 
about  6  h  for  each  case.  When  the  results  satisfied  the  following 
conditions,  the  solutions  were  treated  as  converged: 


in  Table  5.  To  minimize  the  precision  uncertainty  of  a  measurement, 
the  experiment  was  measured  ten  times  in  the  same  way  and  in  the 
same  circumstances,  and  then  the  average  was  taken. 


R  =  \anb^nb  +  b  -  apQ  <  \0  6 

domain 


(21a) 


5.  Results  and  discussion 


Ci+1  -  Ki 


<  io  5, 


(21b) 


where  R  represents  the  residual  sum,  and  £  is  a  general  dependent. 
The  subscripts  i  and  nb  are  the  number  of  iterations  and  the 
neighborhood  grid  points,  respectively. 

In  order  to  verify  the  accuracy  of  the  present  numerical  scheme, 
we  compared  our  1-D  numerical  results  with  that  by  Fraisse  et  al. 
[42]  for  the  1-D  analytical  solution  with  =  500  K  and  Tc  =  300  K, 
as  shown  in  Table  4.  The  discrepancy  is  less  than  2.64%. 


Fig.  6(a)  and  (b)  show  the  voltage-current  ( V—I )  and  power- 
current  (P— /)  curves  depicting  the  comparisons  for  the  TEG  module 
with  and  without  a  plate-fin  heat  sink  (N  =  6,  Hfin  =  30  mm).  The 
experimental  results  are  denoted  by  circular  symbols.  For  the 
convenience  of  verification  of  the  accuracy  of  the  numerical 
simulation,  the  work  fluid  is  treated  as  hot  air  in  the  cases  shown  in 
Fig.  6.  It  can  be  seen  that  the  predicted  numerical  data  for  the  V-I 
and  P—l  curves  are  in  good  agreement  (within  9%)  with  the 
experimental  data.  In  addition,  the  results  show  the  open  circuit 
voltage  Voc  to  be  up  to  2.261  V  (564%  increase)  for  the  TEG  module 


4.  Experimental  setup  and  uncertainty  analysis 


The  schematic  diagram  of  the  experimental  setup  is  illustrated 
in  Fig.  5.  The  wind  tunnel  system  shown  consists  of  a  blower,  12 
annular  heaters  (35  kW),  a  settling  chamber,  the  TEG  modules  and 
an  exhaust  pipe.  The  air  flow  was  driven  by  a  5  HP  centrifugal  fan 
with  an  inverter  to  provide  various  inlet  velocities.  During  the 
isothermal  test,  it  was  determined  that  the  variations  in  the  tem¬ 
perature  in  different  locations  were  within  ±0.2  °C  and  that  a  free 
stream  turbulence  intensity  of  less  than  2%  was  achieved.  The  air 
temperatures  at  the  inlet  and  the  exit  zones  across  the  test  section 
were  measured  using  two  pre-calibrated  RTDs  (pt-100  Q)  that  had 
an  accuracy  of  ±0.05  °C.  The  desired  air  velocity  through  the  test 
section  was  set  by  means  of  adjusting  an  inverter.  The  air  velocity 
was  measured  by  a  hot  wire  with  ±2.0%  accuracy.  Furthermore,  all 
wind  tunnel  systems  were  surrounded  by  insulators  in  order  to 
reduce  heat  loss. 

The  TEG  module  was  cooled  with  water,  and  its  temperature 
was  controlled  at  27  °C  by  a  thermostat  reservoir.  The  water 
volumetric  flow  rate  was  measured  by  a  magnetic  volume  flow 
meter  with  0.002  L/s  resolution.  The  heat  transfer  coefficients  in  the 
water  channel  ( hc )  could  be  estimated  by  using  the  Dittus-Boelter 
equation.  A  high  power  electrical  resistor  array  was  connected  in 
series  to  the  system  which  was  used  to  capture  the  matching  load 
for  the  optimized  power  output.  The  ammeter  and  the  voltmeter 
were  connected  to  measure  the  current  in  the  circuit  and  voltage 
cross  the  external  load  resistor,  respectively.  A  digital  multi-meter 
(Fluke  170)  was  used  to  measure  the  voltage  with  an  accuracy  of 
0.15%  and  the  current  with  an  accuracy  of  1.0%.  All  the  data  signals 
were  collected  and  converted  using  a  data  acquisition  system 
(a  hybrid  recorder). 

When  the  uncertainties  in  measured  and  calculated  quantities 
were  estimated,  the  bias  and  precision  errors  were  both  considered. 
These  elemental  errors  were  combined  to  obtain  an  overall  un¬ 
certainty  in  a  measured  quantity  using  the  root-sum-squares  (RSS) 
method,  which  may  be  estimated  as  follows: 


(22) 


where  is  the  uncertainty  in  the  result  and  wi,  W2,  . . .  ,  wn  are  the 
uncertainties  in  the  independent  variables,  while  W  is  a  given 
function  of  the  independent  variables,  Xi,  X2, ... ,  Xn.  Uncertainties  in 
the  reported  experimental  values  of  the  effective  heat  transfer  co¬ 
efficient  (h)  and  the  power  output  (P)  were  estimated  using  the 
method  suggested  by  Moffat  [43].  The  uncertainties  are  tabulated 


fin  area  (cm2) 


9  50  100  150  200  250  300  350  400 


Hju,  (mm) 

(a)  extra  pressure  drop  (AP)  and  extra  pumping  power  density  (Pfa,/Ab) 
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Fig.  8.  The  AP,  PfanMb,  Pteg  and  r\  vs.  fin  height  (Hfin)  for  different  values  of  flue  gas 
velocity  with  N  =  6  and  Tgas  =  600  K. 


278 


J.-Y.  Jang  et  al.  /  Energy  53  (2013)  270-281 


fin  area  (cm2) 

9  50  100  150  200  250  300  350  400 


Fig.  9.  The  PnetMb  and  PtegM&  vs.  fin  height  (Hfin)  for  different  values  of  flue  gas 
velocity  with  N  =  6  and  Tgas  =  600  K. 

with  a  plate-fin  heat  sink  compared  to  that  of  the  TEG  module 
without  a  plate-fin  heat  sink  (V0c  =  0.401  V).  It  can  be  seen  that  the 
maximum  power,  Pj eg,  occurred  when  the  external  load  resistance 
was  equal  to  the  effective  internal  resistance  of  the  TEG.  A  close 
look  in  Fig.  6(b)  indicates  that  the  Pteg  with  a  plate-fin  heat  sink  can 
be  significantly  increased  up  to  33  times  compared  to  that  without 
a  plate-fin  heat  sink. 

Fig.  7(a)  illustrates  the  temperature  distribution  of  the  first  TEG 
module  with  a  plate-fin  heat  sink  for  three  different  fin  heights 
(Hfin  =  25  mm,  50  mm  and  75  mm)  at  the  number  of  fins  N  =  6  and 
a  flue  gas  velocity  Vin  =  10  m/s.  The  cold  surfaces  are  assumed  to  be 


exposed  to  heat  convection  with  the  heat  transfer  coefficient 
hc  =  1000  W/m2  K,  which  is  decided  with  reference  to  the  tem¬ 
perature  differences  between  an  interior  surface  and  the  water 
temperature  (TWater  =  300  I<).  An  increase  of  fin  height  leads  to 
intensification  of  the  heat  transfer  rate.  The  main  reasons  for  this 
are  that  new  thermal  boundary  layers  are  developed  and  that  the 
heat  transfer  area  is  increased.  In  this  case,  the  predictions  show 
that  the  average  base  temperature  of  the  heat  sink  at  506.4  K, 
527.4  K  and  532.3  I<  occur  at  Hfin  =  25  mm,  50  mm  and  75  mm, 
respectively.  Thus,  a  proper  heat  sink  size  can  indeed  decrease  the 
total  thermal  resistance  due  to  the  presence  of  a  larger  heat  transfer 
area  from  the  hot  gas  side  to  the  cold  water  side.  Additionally, 
Fig.  7(b)and  (c)  show  the  temperature  and  electric  potential  dis¬ 
tributions  of  TEG  for  three  different  fin  heights  (Hfin  =  25  mm, 
50  mm  and  75  mm),  respectively.  It  can  be  seen  that  the  average 
temperature  difference  between  the  top  and  bottom  surfaces  of  the 
TEG  for  Hfin  =  75  mm  is  150.1  K,  which  is  larger  than  the  other  two 
cases  (133.7  K  and  146.9  I<).  Thus,  the  voltage  at  maximum  power 
Vmp  =  T41  V  (11%  and  2.1%  increase)  for  Hfin  =  75  mm  as  compared 
to  the  other  two  modules  Hfin  =  25  mm  and  50  mm,  for  which 
Vmp  =  1.268  V  and  1.383  V,  respectively. 

Fig.  8(a)  and  (b)  present  the  variations  in  the  extra  pressure  drop 
(A P),  the  extra  pumping  power  density  (Pfan/Ab),  the  ideal  power 
output  (Pteg)  and  the  conversion  efficiency  (77)  of  a  TEG  with  the  fin 
height  of  heat  sink  (Hfin),  respectively,  for  three  different  values  of 
flue  gas  velocity  (Vin  =  3,  5  and  10  m/s)  with  the  number  of  fins 
N=  6  and  Tgas  =  600  K.  As  the  plate-fin  heat  sink  was  inserted  into 
the  chimney  tunnel,  although  the  flue  gas  tended  to  bypass  from 
the  clearance  around  the  fin  and  to  redistribute  itself,  it  undoubt¬ 
edly  increased  the  total  tunnel  pressure  drop.  As  expected,  the 
extra  pressure  drop  and  pumping  power  density  were  increased  as 
either  the  fin  height  (Hfin)  or  the  flue  gas  velocity  (Vin)  were 
increased.  As  the  fin  height  of  the  heat  sink  went  from  0  to  100  mm 
for  N  =  6,  the  heat  transfer  area  increased  from  9  cm2  to  387  cm2. 
Thus,  the  ideal  power  output  and  conversion  efficiency  of  the  TEG 
was  relatively  increased.  There  is  an  asymptotic  nature  of  the  Pteg 
vs.  Hfin  curves  when  Hfin  is  greater  than  60  mm.  At  this  time,  the 
conversion  efficiency  was  3.3%,  3.4%,  and  3.5%  for  Vin  =  3,  5  and 
10  m/s,  respectively. 
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Fig.  10.  The  temperature  and  electric  potential  distributions  of  the  first  thermoelectric  module  with  a  plate-fin  heat  sink  for  Hfin  =  30  mm,  Vin  =  10  m/s  and  Tgas  =  600  K. 
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In  order  to  know  whether  or  not  the  net  electric  power  (Pne t)  of 
the  TEG  module  can  actually  generate,  the  extra  pumping  power 
(Pfan)  should  be  subtracted  from  the  ideal  electric  power  (Pteg)  to 
obtain  the  actual  power.  As  shown  in  Fig.  9,  the  effect  of  fin  height 
on  the  net  power  density  (PnetMb)  for  three  different  flue  gas  ve¬ 
locities  ( Vin  =  3, 5  and  10  m/s)  with  the  same  number  of  fins  N  =  6  is 
described.  The  dashed  lines  denote  the  results  for  the  ideal  power 
density  for  a  TEG,  and  the  results  of  net  power  density  are  denoted 
by  solid  lines.  When  the  flue  gas  velocity  is  relatively  slow,  the  net 
power  density  doesn’t  exhibit  notable  differences  with  increases  in 
the  fin  height  of  the  plate-fin  heat  sink.  However,  there  is  a  notable 
difference  between  the  ideal  power  density  and  the  net  power 
density  as  the  flue  gas  velocity  is  increased.  Thus,  an  optimum  fin 
height  Hfin  exists  in  correspondence  with  the  maximum  net  power 
density  PnetAV  As  an  example,  the  optimal  fin  height  Hfin  is  about 
50  mm  for  N  =  6,  Vin  =  10  m/s  and  Tgas  =  600  K. 
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(a)  extra  pressure  drop  (AP)  and  extra  pumping  power  density  (P/atMb) 

fin  area  (cm2) 

65.7  75  100  125  150  175 


N 

(b)  power  output  (Pteg)  and  conversion  efficiency  (i]) 

Fig.  11.  The  AP,  PfanMfa.  Pteg  and  7]  vs.  number  of  fins  (N)  for  different  values  of  flue  gas 
velocity  with  Hfin  =  30  mm  and  Tgas  =  600  K. 


Similarly,  the  performance  of  the  TEG  module  with  a  plate-fin 
heat  sink  for  different  numbers  of  fins  can  be  evaluated.  Fig.  10(a) 
displays  the  temperature  distribution  of  the  first  TEG  module  with 
a  plate-fin  heat  sink  for  three  different  numbers  of  fins  and  a  fixed 
fin  height  Hfin  =  30  mm.  When  the  number  of  fins  is  increased  to 
enlarge  the  heat  transfer  area,  the  heat  transfer  rate  is  relatively 
raised.  It  can  be  observed  that  the  average  base  temperatures  of  the 
heat  sinks  are  486.4  K,  515.4  I<  and  533.5  K  for  N  =  4,  6  and  8, 
respectively.  As  can  be  seen  in  Fig.  10(b)  and  (c),  the  temperature 
and  electric  potential  distributions  of  the  TEG  for  three  different 
numbers  of  fins  (N  =  4,  6  and  8)  are  illustrated  at  a  maximum  TEG 
power.  Respectively,  the  average  temperature  difference  between 
the  top  and  bottom  surfaces  of  the  TEG  are  120.9  K,  138.2  K  and 
150.7  K  for  N  =  4, 6  and  8,  respectively.  The  voltage  values  occurring 
at  a  maximum  power  at  N  =  8  (Vmp  =  1.415  V)  increases  22%  and 
7.4%  compared  to  N  =  4  (Vmp  =  1.156  V)  and  N  =  6  (Vmp  =  1.318  V). 

Fig.  11(a)  and  (b)  demonstrate  the  variations  of  the  extra  AP,  Pfan/ 
Ab,  Pteg  and  r/  with  the  number  of  the  fins  (N),  respectively,  for  three 
different  values  of  flue  gas  velocity  (Vin  =  3,  5  and  10  m/s)  with  the 
fin  height  Hfin  =  30  mm  and  Tgas  =  600  K.  As  the  number  of  fins 
varies  from  4  to  8  in  this  case,  the  heat  transfer  area  increases  from 
84.6  cm2  to  160.2  cm2.  It  is  evident  that  the  extra  pressure  drop  and 
the  extra  pumping  power  density  are  significantly  increased  as  the 
number  of  fins  (N)  is  increased.  Thus,  it  can  be  observed  that  as  the 
number  of  fins  is  elevated  from  4  to  8,  the  value  of  extra  pumping 
power  density  is  increased  up  to  3.6  times  for  Vin  =  10  m/s.  Namely, 
the  extra  pumping  power  increases,  and  thus,  the  power  genera¬ 
tion  costs  are  relatively  raised.  Additionally,  the  ideal  power  output 
(Pteg)  and  conversion  efficiency  (rj)  is  increased  expectedly  as 
either  the  number  of  fins  (N)  or  the  flue  gas  velocity  (Vin)  is 
increased. 

Fig.  12  shows  the  effect  of  the  number  of  fins  (N)  on  PiEG/Ab  and 
Pnet/Ab  for  different  values  of  flue  gas  velocity,  respectively.  As  the 
number  of  fins  increases,  the  net  power  density  increases  as  well  as 
the  ideal  power  density  because  the  increase  in  ideal  power  density 
is  larger  than  the  increase  in  extra  pumping  power  density.  It  is 
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Fig.  12.  The  Pnetl^b  and  Pteg  I ^b  vs.  number  of  fins  (N)  for  different  values  of  flue  gas 
velocity  with  Hfin  =  30  mm  and  Tgas  =  600  K. 
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worth  mentioning  that  the  differences  between  the  net  power 
density  and  the  ideal  power  density  are  enlarged  with  increases  in 
the  number  of  fins.  For  example,  for  Hfin  =  30  mm,  Tgas  =  600  K  and 
Vin  =  10  m/s,  the  difference  between  Pteg/^  and  Pnet/Ab  is  about 
0.011  W/cm2  for  N  =  4,  while  it  is  enlarged  to  0.039  W/cm2  for 
N  =  8.  Furthermore,  as  the  flue  gas  velocity  is  slow  ( Vjn  =  3  and  5  m / 
s),  the  Pnet/Ab  are  very  close  to  the  PrEG/Ab  due  to  the  fact  that  the 
fan  pumping  power  density  Pfan/Ab  is  relatively  smaller  for  the 
lower  flue  gas  velocity  cases. 

Finally,  Fig.  13(a)  presents  the  variations  in  the  ideal  power 
density  (Pteg/^W  of  the  TEG  module  (N  =  6  and  Hfin  =  30  mm)  vs. 


(a)  power  density  of  a  TEG  (PrEG^b) 


Vin  (m/s) 

(b)  net  power  density  ( P„eMb ) 

Fig.  13.  The  PtegM&  and  Pnetl^b  vs.  velocity  (Vin)  for  different  values  of  flue  gas  tem¬ 
perature  with  N  =  6  and  Hfin  =  30  mm. 


the  flue  gas  velocity  (Vin)  for  different  values  of  flue  gas  tempera¬ 
tures  (Tgas  =  500,  600,  700  I<).  The  figure  mainly  shows  the  com¬ 
parisons  with  radiation  effect  and  without  radiation  effect.  It  can  be 
seen  that  the  values  of  PiEG/Ab  with  the  radiation  effect  are  higher 
than  those  without  the  radiation  effect.  As  a  result,  the  P^/Ab  with 
the  radiation  effect  is  improved  by  an  increment  of  7—9%  in  the 
studied  TEG  module  (N  =  6  and  Hfin  =  30  mm)  for  slow  flue  gas 
velocity.  However,  the  increment  is  lowered  with  increases  in  flue 
gas  velocity  because  the  forced  convection  effect  is  gradually 
strengthened.  As  seen  in  Fig.  13(b),  the  comparisons  of  net  power 
density  (Pnet/Ab)  and  the  ideal  power  density  (PrEG/Ab)  vs.  the  flue 
gas  velocity  (Vjn)  for  different  values  of  flue  gas  temperatures  are 
presented.  It  is  evident  that  the  difference  between  PrEG/Ab  and 
Pnet/Ab  is  enlarged  with  increases  in  the  flue  gas  velocity,  and  there 
is  an  appropriate  velocity  for  correspondence  of  the  maximum  net 
power  density.  When  the  flue  gas  velocity  (Vin)  is  increased  over 
6  m/s  at  Tgas  =  700  K,  there  is  greater  pumping  power  penalty  over 
the  raised  power  density,  and  therefore,  the  net  power  density  may 
be  decreased. 

6.  Conclusions 

The  characteristics  of  heat  transfer  behavior  and  power  gener¬ 
ation  of  a  TEG  module  with  a  plate-fin  heat  sink  are  investigated 
numerically  and  experimentally,  and  the  radiation  effect  of  flue  gas 
is  also  considered  in  this  study.  The  maximum  power  Pjeg  of  a  TEG 
module  with  a  heat  sink  can  be  increased  up  to  several  dozen  times 
as  compared  to  a  TEG  module  without  a  heat  sink.  The  numerical 
results  for  the  power  vs.  current  (P— /)  curve  are  in  good  agreement 
with  the  experimental  data  within  an  error  of  9%. 

The  results  indicate  that  the  extra  pressure  drop  and  pumping 
power  density  are  increased  with  increases  in  fin  height  (Hfin)  and 
the  number  of  fins  (N).  Similarly,  both  the  ideal  electric  power 
(Pteg)  and  conversion  efficiency  (rj)  are  also  relatively  increased 
with  increases  in  fin  height  (Hfin)  and  the  number  of  fins  (N).  It  is 
noteworthy  that  the  optimal  fin  height  Hfin  is  about  50  mm  for  the 
correspondence  of  maximum  net  power  density  (Pnet/Ab)  at  N  =  6, 
Vin  =  10  m/s  and  that  Tgas  =  600  K. 

Additionally,  the  PrEG/Ab  with  radiation  effect  increases  7-9%  as 
compared  to  the  PrEG/Ab  without  radiation  effect  in  the  studied  TEG 
module  (N  =  6  and  Hfin  =  30  mm)  for  slow  gas  velocity.  However, 
the  difference  between  the  two  is  smaller  with  increasing  flue  gas 
velocity  because  the  forced  convection  effect  is  gradually 
strengthened.  The  proper  flue  gas  velocity  for  correspondence  of 
maximum  Pnet/Ab  exists  due  to  the  greater  pumping  power  penalty 
as  compared  to  the  raised  power  output.  The  results  obtained  in 
this  study  should  provide  useful  information  for  an  industrial  TEG 
module  design  for  waste  heat  recovery. 
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